arXiv:1508.07972v2 [cond-mat.mes-hall] 19 Oct 2016 


Dissipation and resonance frequency shift of a resonator magnetically coupled to a 

semiclassical spin 

J. M. de Voogd,* J. J. T. Wagenaar, and T. H. Oosterkamp 
Kamerlingh Onnes Laboratory, Leiden University, 

PO Box 9504 , 2300 RA Leiden, The Netherlands 
(Dated: October 21, 2016) 

We calculate the change of the properties of a resonator, when coupled to a semiclassical spin 
by means of the magnetic field. Starting with the Lagrangian of the complete system, we provide 
an analytical expression for the linear response function for the motion of the resonator, thereby 
considering the influence of the resonator on the spin and vice versa. This analysis shows that the 
resonance frequency and effective dissipation factor can change significantly due to the relaxation 
times of the spin. We first derive this for a system consisting of a spin and mechanical resonator and 
thereafter apply the same calculations to an electromagnetic resonator. Moreover, the applicability 
of the method is generalized to a resonator coupled to general two and more level systems, providing 
a key to understand some of the problems of two level systems in quantum devices. 


Resonators and spins are ubiquitous in physics, 
especially in quantum technology, where they can be 
considered as the basic building blocks, as they can 
collect, store and process energy and informationb^. 
The validity of this information is, however, of limited 
duration as these building blocks leak practically always 
to the environment, which on its own can be seen as 
a bath of resonators and spins'”'. If in particular we 
focus on the situation where a resonator is coupled 
to a certain spin, then the spin’s interaction with 
the environment naturally causes, besides a shift of 
resonance frequency, an extra dissipation channel for the 
resonator. Despite this simple qualitative explanation 
and many experimental'^^^ and theoretical efforts®^'^^, 
an applicable full picture that quantitatively describes 
the response of a resonator coupled to a spin and their 
environments is still lacking. Here we derive classically 
the linear response function of the non-conservative 
system consisting of a resonator and a semiclassical spin. 
We show that the quality factor and resonance frequency 
of the resonator can be significantly influenced due to 
the relaxation times of the spin. 

We start with a Lagrangian description, that includes 
the degrees of freedom of the resonator and the spin, to 
find the coupled equations of motion (EOMs) that de¬ 
scribe the resonator displacement and the spin magnetic 
moment, finding that this magnetic moment depends 
on the path the resonator takes. This is fundamentally 
different from conventional magnetic force microscopy 
(MFM)'^^, where one assumes a fixed polarization of the 
spins, like in magnetized samples. Even in magnetic 
resonance force microscopy (MREM), which is usually 
focused on paramagnetic spins, it is generally assumed 
that the spin is not, or at least not significantly, influ¬ 
enced by the resonator'^®^^'®. We will show that this 
influence actually opens the dissipation channel and 
that the resonance frequency shift is more subtle than 


* voogd@physics.leidenuniv.nl 



FIG. 1. Schematic representation of spin fj, interacting with 
two types of resonators, a) The mechanical resonator with 
spring constant k and displacement q{t) of the magnet. The 
dashed line shows the position axis that is used in figure 2. 
b) The electromagnetic resonator as a lumped element device. 
The current I (t) through the inductor L changes the magnetic 
field at the position of spin fi. 


generally assumed. 

Eurthermore, we find in our analytical results that the 
interaction amplitude as function of temperature is a 
curve that for certain conditions shows an optimum, 
see Eig. 3, similar to the curves found in experiments 
where the tails have heuristically been fitted with power 
laws'®’®. Parts of the analysis we present here have been 
used by den Haan et al .to explain the experimental 
results obtained by approaching a native oxide layer 
on silicon with an ultra-sensitive MREM probe. The 
equations derived in this paper were found to closely 
resemble the measured shift in resonance frequency and 
reduced quality factor as function of temperature and 
resonator - spin surface distance. 

Although we start calculating the susceptibility of 
the more intuitive mechanical resonator, we will as 
well derive explicitly the (additional) impedance for 
electromagnetic resonators (see Fig. lb versus la, and 
Sec. IV), thereby making the results suitable for direct 
use by other fields in physics. Moreover, we will show the 
applicability of the theory to the case of the resonator 
coupled to two level systems (2LSs) and higher level 
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quantum systems. 

Finding an accurate description of the interaction 
of the building blocks of quantum devices with the 
environment can be seen as a widespread and major 
research area since not being able to understand, control 
and minimize the interaction is a major bottleneck in: 
the field of quantum computing^d i ^ detector fabrication 
in astronomy^^ds^ MRFM and high resolution 
and the development of optomechanical-like hybrid 

quantum devices" ^’^ 2 ^ 


I. BASIC PRINCIPLES 

The configuration of our theoretical analysis is given in 
Fig. la. A semiclassical spin, with magnetic moment fi, 
is located at laboratory position rg and feels a magnetic 
field B{rs, t) that is produced by a magnet. The magnet 
is attached to a mechanical resonator that has spring 
constant k and (effective) mass m. The origin of the 
laboratory frame is chosen to be the equilibrium position 
of the magnet’s center. The displacement of the magnet 
from this equilibrium position is denoted by q(t). See 
Fig. 1. The Lagrangian for this system is given by 

L = ]^mq^-]^kq^ + B{q) + Is. (1) 

Is stands for an expression with the internal spin degrees 
of freedom that needs to be included to derive the spin 
EOM. A more detailed account is left in App. A. 

The resonator-spin system does not live in an isolated 
world. Therefore we include dissipation and decay to the 
environment into the EOMs. The hrst differential equa¬ 
tion, derived with respect to the resonator displacement, 
includes the Raleigh dissipation —yq of the resonator. 
This results in 

d 

mq + ^q + kq-^i-—B = Fext{t), (2) 

oq 

where the last term, Fextit), is an external force that is 
exerted on the resonator. 

Starting with the Lagrangian, which contains the degrees 
of freedom for the resonator and the spin, leads to the 
force interaction term —/x • ^B. This is the same 
as —/X • Vi?||q, because of the vanishing curl of the 
magnetic field in free space. Here is the gradient 

of the magnetic held component in the direction of the 
movement of the resonator. In MRFM —/x-is often 
derived from calculating the force-held from the gradient 
of the potential energy V (/x • B), assuming that /x does 
not depend on the position of the resonator^^. However, 
as follows the classical path, we will show by solving 
the spin EOM that fj, is inhuenced by the resonator and 
it is therefore a priori not at all obvious that = 0 as 
long as the spin degrees of freedom are not dehned. 


The other set of differential equations can be found 
by deriving the EOM with respect to the spin degrees 
of freedom. Since the spin interacts with the environ¬ 
ment, we can expect an effectively decaying amplitude 
that is often described by Ti and T 2 ; the time constants 
associated with the decay of the semiclassical magnetic 
moment longitudinal and perpendicular to the magnetic 
field, respectively^"^. If one assumes that the system con¬ 
sists of an ensemble of paramagnetic spins, instead of 
one,the average magnetic moment per spin decays to a 
certain equilibrium vector /Xqoj according to the master 
equation"’’. However, if a single spin over time has on 
average the same behavior as the average of an ensem¬ 
ble at a certain moment, i.e. the spin satisfies ergodicity, 
then we can combine the ensemble’s master equation and 
the single spin EOM to hnd a differential equation that 
describes the average behavior of the single semiclassical 
spin. This is the Bloch equation: 

lil = 'YsfJ. X B +T~^ {flao - II) ■ (3) 

Here 7 s is the gyromagnetic ratio and T~^ = 

(^1 — BB'^'^ + t^BB^, where the hat denotes the 

unit-vector in the direction of the specified vector. 

The spin equilibrium magnetic moment fiao (t) is the vec¬ 
tor to which the spin magnetic moment would decay to 
if given the time. As the resonator moves, the magnetic 
field changes, and so does We will assume that the 
environment of the spin is a heat bath, connected to the 
spin by means of the relaxation times. However, does 
the spin’s equivalent spin ensemble have a well defined 
temperature? As derived in the original paper of Bloem- 
bergen et al. , the differential equation describing the 
population difference n for particles in a two level system 
is 


dn 

dt 


-2Wn + 


uq — n 
Ti 


( 4 ) 


where W is the probability rate that the particle changes 
energy level due to an applied field and no is the popu¬ 
lation difference between the energy levels when the en¬ 
semble has the temperature of the heat bath. In other 
words —2ITn is proportional to the incoming energy and 
is the connection to the heat bath. This results in 


I + 2WT1’ ’ 

where n^o is the steady state solution. Thus when 
2WTi ^ 1 the spin ensemble, and hence our semiclas¬ 
sical spin, is connected well enough to the heat bath to 
assume that our spin has a well defined temperature. For 
spin-4 this condition yields^® 

9"irig(a;)«I, (6) 


where B' = 



and g (w) the spin’s normalized 


absorption line that is usually described by a Lorentzian 
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FIG. 2. This graph shows the single spin contribution to the 
spring constant as function of a position axis parallel to the 
direction of resonator movement, as visualized by the dashed 
line in figure la). In the simulation we attached a magnetic 
dipole (with magnetic moment of 19 pAm^ in the direction of 
q) on a mechanical resonator. The resonator is connected, by 
means of the magnetic field, to an electron spin at a tempera¬ 
ture of 300 mK. The distance between the center of the dipole 
and a; = 0 is 2.5 pm. To demonstrate the spatial behavior of 
the K-terms we avoided imaginary terms by setting Ti = 0 in 
K 2 and T 2 = 0 in K 3 . The solid line shows the sum of these 
K-terms. 



or Gaussian that peaks around the Larmor frequency. 
This makes this condition hard to satisfy when the res¬ 
onator has a resonance frequency around the Larmor fre¬ 
quency, and one should minimize the resonator’s move¬ 
ment q. When this condition is not met, the spin sat¬ 
urates and the temperature increases or might be un¬ 
defined^®. However, for example in MRFM, mechanical 
resonators tend to have resonance frequencies much lower 
than the Larmor frequency and so it is much easier to sat¬ 
isfy this condition. 

Assuming the condition is satisfied we can now derive 
fXao from the canonical ensemble and find for spin-^ 

Moo =/istanh(,5ps (7) 

where /? = is the inverse temperature and = Shjs 
is the magnitude of the non-averaged spin magnetic 
moment with spin number S = ^. This result can easily 
be generalized for other spin numbers as is done in 
App. B. For simplicity we will stick to the formula for 
spin-i particles here. 


term up to linear order in q. Higher order terms will 
give rise to nonlinear effects. Interaction terms with 
even powers in q are usually experimentally uninter¬ 
esting since they will produce even multiples of the 
fundamental resonance frequency. These multiples are 
not measured or can easily be filtered. Uneven powers 
of q can, however, lead to disturbing nonlinear effects 
like Duffing^^. One can lower the amplitude of q to 
suppress higher order terms and therefore the nonlinear 
effects, but in experiments this is usually limited by the 
signal-to-noise ratio. 

The zeroth order term does not contribute to the 
dynamics of the system, however it does give rise to a 
constant deflection of the resonator. This can be solved 
by shifting the origin of the laboratory frame by the 
amount of the deflection; this causes, however, a (usually 
small) change of the coordinates of the spin. We will 
provide an estimate of the deflection in App. C and 
leave it further out of account. 

To find the interaction term —/jl ■ ^ up to first 
order in q, we need to solve Eq. 3 and find the constant 
and g-dependent parts. By substituting g —>■ Ag we use 
perturbation theory to find 



Mo • + A (Ml • B' - g/xo • B") + O (A") , 

( 8 ) 


where B' = 




was defined previously and B" = 


Here m is perturbed into a g-independent 


part Mo Eind a linear term fii. The higher order terms 
O (A^) can be omitted, as well as the first term on the 
right hand side that only gives rise to the constant de¬ 
flection. 

At first we are mostly interested in solutions that do 
not decay over time and do not depend on initial condi¬ 
tions because then the linear response function can con¬ 
veniently be given in the Fourier domain which makes it 
easy to compare with experiments. The Fourier Trans¬ 
form T{ } of the linear response function, or simply sus¬ 
ceptibility X (w) = , can be calculated from Eq. 2 


x(w) 


I 

k — muj"^ -b fqw -I- k, ’ 


( 9 ) 


II. SUSCEPTIBILITY 


To find the resonance frequency and quality factor of 
the resonator, we will need to calculate the interaction 


where k = K 1 + K 2 + K 3 , with ki = -fig -B" and K 2 + K 3 = 
^ • Appendices C,D present the calculation of the 















4 


K-terms, which turn out to be: 


Ki = 

K2 (w) = 
K 3 (w) = 


tanh (/3/rsBo) 
l3/j,sBo 


Bq cosh^ (/3/is-Bo) 

--^tanh i/3fj.sBo) 
Bq 


B' 




1 


1 + iujTi ’ 


1 - 


2a-(a;r2f+ia;r2(l + 2a) 
(1 + iujT2f + {iO,T2f 


( 10 ) 

( 11 ) 


( 12 ) 


where Bq = B {q — 0) and the notation and 
is used to indicate the part of v parallel and perpendic¬ 
ular to Bq respectively for any vector v. K 2 and kq are 
derived from rti,, - and Ui - respectively. 

II-Sq -L-Bq 

If we compare this result with the conventional approach 
that neglects the effect of the resonator on the spin, we 
see that in that approach we have only the term 
However, ki is real and therefore it cannot describe the 
extra dissipation channel that has been seen in experi¬ 
ments^®. The derivation which has been done here does 
include the linear effect of the resonator on the spin and 
vice versa. This produces two extra terms in the linear 
response function that are partly imaginary. Each of the 
K-terms is shown separately in figure 2 as a function of 
the spin position. This position axis is indicated in fig¬ 
ure 1 by the dashed line. Which effect these terms have 
in practice, where usually more than one spin is present, 
will be shown in the next section. 


III. SPIN BATH - RESONATOR COUPLING 


We assume that all spins in the system act individ¬ 
ually and do not influence each other, except through 
the relaxation times. We can then sum over the k- 
terms for each spin to find the susceptibility of the 
resonator connected to a whole ensemble of spins, i.e. 
K = J2s ^i(’'s) + «^ 2 (»’s) + K 3 {rs). Moreover, if the spins 
in the sample have an average nearest neighbor distance 
smaller than the typical spatial scale of the applied mag¬ 
netic field, we can see the sample as a spin continuum 
and hence, instead of summing, integrate over the sam¬ 
ple with spin density p{r). 

If we calculate the result for a volume with constant spin 
density, it can be found by partial integration of the vol¬ 
ume in the direction of the movement of the resonator 


k(w) =pf3plC 


{ioTif + iujTi 
l + (c^Tl)2 


boundary term -|- O 



1 

-w2)r| 


(13) 


with 


C= / 


B!r 


Iv cosh^ {PpsBq) 


a) 



FIG. 3. Calculated frequency shift and added dissipation of 
a mechanical resonator dne to dangling bonds on a silicon 
surface, equivalent to the setup of den Haan et al. . a) Im¬ 
pression of a NdFeB magnet (with magnetic moment 19 pAm^ 
in the direction of q) attached to an ultrasoft silicon cantilever 
with spring constant fc = 70 pN/m, together leading to a nat¬ 
ural frequency of = 3 kHz. The center of the magnet is 
positioned at a distance of 2.2 pm to the silicon sample. The 
surface of the sample has a native oxide containing 0.14 elec¬ 
tron spins/nm^ that are visualized by the red balls (not to 
scale). The graphs b) and c) show the resonance frequency 
shift and the damping of the cantilever. The results are shown 
for various Ti, showing a maximal opening of the additional 
dissipation channel for Ti = l/a;o. 


The boundary term vanishes when the volume bound¬ 
aries in the q-direction are large. The O 


can be neglected for resonance frequencies away from the 
Larmor frequency and for r 2 ^ . 

From K we can calculate the frequency and Q-factor 
shifts as seen in experiments by den Haan et al. For 
Qq = ^ ^ the susceptibility has a maximum 


around the natural frequency ujq = y Then, as long 
as the influence of the spin leads only to a small correction 
of the susceptibility, i.e. k ^ k, the relative frequency 
shift is given by 


(14) 


Aw ^ 1 Re (k(wo)) 

ujq 2 k 


(15) 
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T(K) 


FIG. 4. Simulation of frequency shift and added dissipation 
of an electromagnetic resonator due to dangling bonds at the 
sample’s surface, a) Impression of an RLC-circuit with 10 
GHz natural frequency and 0.25 nH inductance that consists 
of a 50 qm x 50 qm square which is positioned 50 nm above 
a surface with 0.14 electron spins/nm^ b,c) Galculated results 
for a static external magnetic field of 0.1 T that is oriented 
out of plane (solid curve) and in plane (dashed curve). For 
this simulation we assumed T 2 = 0.01 qs. 


The imaginary part of k causes the change in Q-factor. 
The new Q-factor is given by 

11 Im(K(a;o)) 

o'oo’'— k —■ 

In Fig. 3 we show an example of an experiment with a 
magnet attached to an ultrasoft cantilever, which is posi¬ 
tioned above a silicon sample. The native oxide contains 
electron spins that interact with the resonating magnet. 
The frequency shift and quality factor depend differently 
on Ti. In this simulation we have set T 2 to zero only af¬ 
ter we checked that the O term in Eq. 13 can indeed be 
neglected: setting T 2 = Ti gives an additional frequency 
shift of about 1 nHz and a five orders of magnitude lower 
shift in Q-factor compared to the results shown in Fig. 3c. 


IV. SPIN - ELECTROMAGNETIC RESONATOR 

In this section we calculate the complex impedance 
coming from a spin interacting with an electromagnetic 
resonator. The derivation is very similar to the mechan¬ 
ical resonator and hence we will largely copy the results. 
We will assume that the system can be described by a 


lumped element model, which is a valid approximation 
when the typical size of the system is much smaller than 
the wavelength. The results might be generalized to work 
for other resonators by using the distributed element 
rnodek’^^^. However, this can become rather complicated 
depending on if it is necessary to calculate the interac¬ 
tion between resonator and spin using the retarded time 
(Jefimenko’s equations). Moreover, it could be that the 
interaction depends on the current density rather than 
the current, all of which is outside the scope of this pa¬ 
per. We conveniently describe a series RFC circuit, see 
Fig. lb. 

As there is a direct analogy with the mechanical res¬ 
onator, it is straightforward to write down the complete 
Lagrangian and derive the EOM. From this we calculate 
something similar to the susceptibility, but more com¬ 
monly used in electromagnetism, the impedance Z{ijj) = 

AM 

The electromagnetic analog of the displacement q is the 
charge Qe- However, instead of writing down Qe and 
‘momentum variable’ Qe, we prefer to work with the cur¬ 
rent I = Qe- The ‘position variable’ Qe then becomes 
JdtI. This results in the RLC-resonator’s EOM as 

+ (17) 


The resulting interaction term is slightly different com¬ 
pared to that of the mechanical resonator. The zeroth 
order term vanishes conveniently due to the time deriva¬ 
tive, leading to the impedance interaction term z(uj) = 

The spin’s EOM does not change, apart 
from change of variable q ^ I. This results in an extra 
impedance z = zi + Z 2 + z^, equivalent to the K-terms, 
where 


zi = iujfis tanh [P^sBq) 


B 


Z 2 = 


Bo cosh^ (PusBq) 


.Bo 

B' 


1 

1 -I- iujTi ' 


-^3 = tanh {PfisBo) 

Bq 


B'. 


±Bo 


2^ - (ujT2f + iujT2 (1 + 2^y 
(1 + iujT2f + {uJsT2f 


(18) 

(19) 


( 20 ) 


The resonators complex impedance then becomes 

Z{uj\ = iuL H. -\—^—— -\- z. 

iujC 


( 21 ) 


It is much harder to simplify the z-terms as done in 
Sect. HI when partially integrating over a whole sample 
because I is, unlike q, not a Cartesian direction. How¬ 
ever, one thing simplifies the z term reasonably: the law 
of Biot-Savart shows a linear dependence on I implying 
that Zi vanishes. Note that it is very well possible that 
the frequencies of interest are comparable to or Wg. In 
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this case one should calculate the whole term. Moreover 
one should be careful with the implied condition of Eq. 6, 
i.e. /2ri5( w) <C 1 when probing the resonator. 

In Fig. 4 we provide an example of an electromagnetic 
RLC-circuit fabricated on top of a silicon sample with a 
native oxide. The electron spins inside the native oxide 
couple to the inductor changing the resonators resonance 
frequency and Q-factor. 


V. RESONATOR COUPLING TO OTHER 
SYSTEMS 

So far we have done nothing more than rigorous math 
to calculate the susceptibility of a system were the phys¬ 
ical process is precisely known. However, the physical 
nature of the interaction between a resonator and a gen¬ 
eral two level systems (2LSs) can be different from the 
simple magnetic field interaction and will often even be 
unknown. This subject has been studied in glassy sys¬ 
tems long before it found its application in quantum tech¬ 
nology 'h The field revived when it was found in exper¬ 
iments that the electric permittivity and loss factor of a 
nonmagnetic glass do actually depend on the magnetic 
field'^^. It was only until recently, around the same time 
as this paper appeared on a preprint server, that Jug 
et al. provided an intuitive and elegant explanation 
based on a JB-field dependent density of states and heat 
capacity. Indeed, expanding the average energy term, as 
we did in App. B, leads to the heat capacity which re¬ 
sulted in K 2 oc ^2 oc with x the Zeeman energy — 

temperature ratio. These similar results in combination 
with the results obtained in this paper imply two things: 
First the B does not have to be the physical magnetic 
field. It is always possible to rewrite the two state Hamil¬ 
tonian to 

H = Eo + ^(T-{Ba-qB' + q^B" + ( 22 ) 

where B can be any field that splits the energy levels, 
leading to an energy difference e when q = 0. Here Eq is 
an uninteresting energy-offset and cr is a vector contain¬ 
ing the Pauli matrices. The interaction strength is deter¬ 
mined hy -^B, hence it is important that B depend on 
q, which is the generalized coordinate of the mechanical 
resonator, or generalized velocity of the electromagnetic 
resonator. Because the expectation values of the Pauli 
matrices cr are described by the Bloch equations, the 
derivations in this paper apply to any resonator-2LS sys¬ 
tem. Just substitute /j.^ —into the k and z terms. 
Secondly, this result can be easily generalized to a system 
with 2S + I energy levels (with S an integer or half in¬ 
teger) by expanding the Brillouin function from App. B 


and substituting 

tanh {P^sBq) —)• {2S + I) coth ((S' -|- |)/3e) — coth 

(23) 

PfJ-sBo -^(2S -b l)^/3e \l3e 

cosh^ (/J^gHo) sinh^ ((S -I- |)/3e) sinh^ (i,0e) 

(24) 

into the k and z-terms. This 2S -I- I-state quantum 
system must be isomorphic to a spin-S particle and 
hence meet two conditions: I) the energy levels are 
equally spaced and 2) transitions are only possible to 
adjacent energy levels. 


VI. DISCUSSION AND CONCLUSIONS 

We have calculated the linear response function of a 
mechanical and electromagnetic resonator coupled to a 
spin. The linear response function of the resonator shows 
extra terms that result in a shift of the resonance fre¬ 
quency and a drop of the Q-factor of the resonator, com¬ 
pared to the bare resonator characteristics. Moreover, 
we have generalized these results to the coupling with an 
energy level system with an arbitrary amount of equally 
spaced energy levels. In practice this means that despite 
having nonmagnetic samples and frequencies that are not 
even close to the Larmor frequency, one encounters dissi¬ 
pation of the resonator due to the inhomogeneous field it 
creates. Eventually this might not be a surprise since the 
resonator alters the heat capacity of the spin’s equivalent 
spin ensemble. Although this is closely related to the 
magnetic loss enhancement in nonmagnetic glassy sys¬ 
tems"*'"*, we did not find any description in literature that 
provides a quantitative and detailed account of how this 
influences the linear response of the resonator, despite 
the many reported and unexplained results’'^” ‘ . The re¬ 
sults presented here have been experimentally verified*® 
and have been used to calculate the frequency shift in a 
simple, yet powerful, saturation measurement protocol'^'*. 
We have chosen to do the calculations completely in the 
(semi)classical regime as we are especially interested in 
the expectation value of spin and resonator. Moreover 
this leads to an intuitive description and fairly simple 
calculations. The classical treatment has it limitations 
though: Berman et al. have raised the point that, if 
the cantilever position is constantly measured, there is an 
influence on the spin because of the projections that are 
constantly occurring in the act of measuring. This might 
introduce random quantum jumps which, when they are 
not time averaged over timescales longer than Ti, are not 
taken into account in our description. Furthermore, when 
pulses are applied, for example in spin resonance tech¬ 
niques, a precise time evolution of the system is needed. 
Moreover, sending hard pulses might violate the condi¬ 
tion for the temperature and linear response of the spin 
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that we have encountered in Sec. I. In this case one might 
move to a calculation involving the spin-operators. The 
theory presented here would still give a fair indication 
about the enhancement of dissipation, which is of impor¬ 
tance in the field of hybrid quantum systems that are 
pushing the limit of macroscopic superpositions^^’^^. 


single phonon nanomechanics project of the Foundation 
for Fundamental Research on Matter (FOM), which is 
part of the Netherlands Organisation for Scientific Re¬ 
search (NWO). 
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Appendix A: Resonator - Semiclassical Spin Lagrangian 

The semiclassical magnetic moment can be seen as a vector with an azimuth (j) and a polar angle 9, where the 
poles of the spherical coordinate system {0 = 0° and 180°) lie on the axis parallel to the magnetic field. 9 and (j) can 
be seen as the two degrees of freedom that a spin has. Then the Lagrangian L = ■ B{q) + Sh(f> cos 9 reveals the 

Bloch equations for a spin-S' particle, but then without decay and for magnetic moment instead of magnetization. 
The last term of the Lagrangian describes the internal dynamics of the spin. Substituting this into the full Lagrangian 
in Eq. 1, we find 

L = — ^kq^ + H ' B{q) + Sfi^cos9. (Al) 


Appendix B: Equilibrium magnetic moment 

By definition of the equilibrium vector we can state that —Roo • B = (E), where (E) is the equivalent ensemble 
average for the energy, or for a single spin the averaged energy over all the points in time with equal q. The limited 
energy levels make it easy to calculate the average energy: For spin-S there are 25" J- 1 energy levels with energies 
Ek = —kgsfJ-s |.B| with k = —S, —S J- 1,..., S'. Using the relation between internal energy and the canonical partition 
function, this results in 

Moo = (2S -I- l)coth ((2S -I-1) /3fis \B\) - coth (/3/is \B\)'^B (Bl) 

g,si&nh{f3fj,s\B\) B. (B2) 

This result is also known as the Brillouin function for the Zeeman energy. The imposed direction B follows from 
Curie’s law. The result might be different when the spin has (strong) interaction with its neighbors and when this 
leads to anisotropic effects, although some of these effects might be included in the q independent part of B. 


Appendix C: Zeroth order solution 

If the magnetic field generated by the oscillating magnet is given by B{r) in the magnet’s rest frame, then in the 
laboratory frame the magnetic field is B{r — Xq). Around the spin position Vg the magnetic field is 

B = Bo-qB' + ^q^B" + .... 


(Cl) 
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Here Bo = B{r,), B'= ^ 


and B” = ^ 2 ^ 


Next we substitute q ^ Xq and expand /Xoo for spin-f up to first order in A and omit higher order terms 


Me 


= Us tanh (PfisBo) Bq - q i tanh (/3/isHo) P± 


cosh^ (PfisBo) 


B' 


(C2) 


Pi I = Bo Bn and 


where Py and P_l are projections parallel and perpendicular to the Bq field respectively, i.e. -ry = -D 0-05 
Pj_ = 1 — BqBq . We also set g —>■ Ag into Eqs. 3 and Cl and set A 0 to get the differential equation to solve for 
fJ-Q- 


Ao = (^IsBox - 7^P± - if i^f^sBo) Bo, 


(C3) 


where the x subscript denotes an antisymmetric matrix such that A^v = v x A for any vector v and A. 

Let M{s) = e“®‘/i.(t) dt be the Laplace transform of the magnetic moment and apply the necessary linear 

algebra to get 


Mo(s) = 


S + J Pj_ + LUsB, 


'Ox 


A fis 


(® + 


s + 7 ^ I \s Ti 


tanh iPusBo) Bo + /x( 0 ) ) , 


(C4) 


with ujs = "fsBo- The inverse Laplace transform yields the general solution for /xq in the time-domain: 

(e~*A2 cos{uJst) —e~*A2 sin(a;st) 0 \ 


/xo(t) =/Xg ("l — e tanh (/3/x,jPo)-B q + I e g cos(a;st) 0 ] /x(0). (C5) 

^ ^ '0 0 e-^/B j 


0 


To retrieve some intuition for the results we choose to present the last term as a matrix which is given in a non-rotating 
Cartesian basis with z = Bq. 

To estimate the static displacement we use /xo(oo), which is of course the same as /Xoo (9 = 0), to find the change of 
equilibrium position 

^ ^ WiJ'sBo) Bo ■ B', (C 6 ) 

where in the last step we neglected 5k, the effective extra stiffness coming from the terms linear in q. 


Appendix D: First order solution 


As argued in the main text, we can ignore the terms that decay or depend on initial conditions. As a consequence we 
can take /xq = ^s tanh (/^/x^Po) Bo- This leads immediately to one of the interaction terms. Taking P {—g/xo ■ B"} = 
Kiqiuj) with q(uj) = P{q{t)} we arrive at 


Kl = -fls 


B': 


tanh i^fisBo ), 


(Dl) 


where 


B': 


= B” ■ Bo 


Next, we need to find /xi. 
only. We find 


Again this is done by substituting q ^ Xq and extracting the terms that are linear in A 


Ai 


(^7s^0x - ^P± - ah 

ftanh(/3^sSo)Pj_-t 
Bo-Li \ 


Mo 


PUsBq a 
cosh^ {Pfj.sBo)) 


B', 


(D2) 
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where C = (^BoB''^Pjl + P±B'B^'^. 

The first line is the same as in Eq. C3 and therefore leads to the same matrix as in Eq. C4 using the same non-rotating 
Cartesian basis with 2 : = Bq. This leads to 


Mi(s) 



P± + ^sBox 



+ w 


2 

S 


+ 



(^tanh{(3^lsBo) 




P±B' - lOsBo X B' 


PfJ^sBp 

cosh^ iPusBo) 



fJ-s 

BqTi 


Q{s), 


(D3) 


with -Mi(s) and Q{s) being the Laplace transform of Hi{t) and q{t) respectively. 

Ml, and thus fii, can be easily split in a part that is parallel and perpendicular to Bp. It follows from Eq. 8 that 
we need specifically the product /xi • B' for the interaction term. So let us write P {/xi • B'} = q(uj) (k 2 + Kp) where 
K 2 and Kp come from the parallel and perpendicular parts of fii respectively. Einally we move to the Fourier domain, 
which is possible since all poles lie in the Re(s) < 0 regime. This leads to 


K2 = - 


Bp 


B', 


fi^'sBp 


cosh^ {f3fj,sBp) 1 + iioTi ’ 


where 


B' 


II Bo 

For Up we find 


= B'^^PiB' 


fJ-s 

Kp = - — 

Pq 




tanh{/3nsBp) 


1 - 


2^-iLoT2f+iLoT2(l + 2^y 
(1 -|- iUjT2) + {uJsT2) 


(D4) 


(D5) 


where 


B' 


J-Bo 


B'^P^B'. 
















